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Using ab-initio calculations we predict the existence of one-dimensional (ID), atomically confined 
plasmons at the edges of a zigzag M0S2 nanoribbon. The strongest plasmon originates from a 
metallic edge state localized on the sulfur dimers decorating the Mo-edge of the ribbon. A detailed 
analysis of the dielectric function reveals that the observed deviations from the ideal ID plasmon 
behavior result from single-particle transitions between the metallic edge state and the valence- and 
conduction bands of the M0S2 sheet. The Mo- and S-edges of the ribbon are clearly distinguishable 
in calculated spatially resolved electron energy loss spectrum owing to the different plasmonic prop¬ 
erties of the two edges. The edge plasmons could potentially be utilized for tuning the photocatalytic 
activity of M0S2 nanoparticles. 


Plasmons are collective electronic excitations that cou¬ 
ple strongly to external fields such as light or fast prop¬ 
agating charges. In metals, their excitation gives rise to 
charge oscillations which can be localized either in the 
interior of the metal a or at its surface. Localized sur¬ 
face plasmon resonances (LSPR) in metal nanoparticles 
have been utilized for chemical sensing jj], cancer treat¬ 
ment [2] and to increase solar cell performance j3j. Two- 
dimensional (2D) plasmons are found at metal surfaces 
and thin films [4], in graphene nn, and other atomi¬ 
cally thin 2D materials such as metallic transition metal 
dichalcogenides |9j. 

Compared to their 2D counterparts one-dimensional 
plasmons have been much less studied. So far, true atom¬ 
ically confined ID plasmons have only been observed 
in self-assembled atomic chains on semi-conducting sub¬ 
strates mm One peculiar fundamental aspect of ID 
metals is the break down of Fermi-liquid theory and the 
transition to Luttinger liquid behavior nans]: At low 
energies the spectral weight of the quasiparticle peak ap¬ 
proaches zero and the life-time broadening becomes of 
the same order as the excitation energy p~4l IT5] . This 
implies that for (ideal) ID metals, the plasmons adopt a 
special status as the only excitations which can couple to 
external fields. 

Recently, the interest in novel 2D materials, such as 
graphene and monolayer M 0 S 2 , has been accelerated by 
the prospects of utilizing their unique electronic prop¬ 
erties for nanoscale (opto) electronics [16, [17]. While 
graphene has shown great potential as a material for 
terahertz plasmonics [18], the finite band gap of M 0 S 2 
makes it unsuitable for applications within this field. In a 
rather different context, M 0 S 2 nanoparticles are used in¬ 
dustrially as hydrodesulphurization catalysts [IT and are 
considered as promising alternative to Pt as low-cost cat¬ 
alysts for hydrogen evolution ED] . The high catalytic ac¬ 
tivity of the M 0 S 2 nanoparticles has been directly linked 
to the presence of metallic edges states on the other¬ 
wise semi-conducting M 0 S 2 nanoparticles [21 . Very re¬ 


cently, an observed extraordinary photo-catalytic behav¬ 
ior of metallic nanoparticles has been shown to originate 
from the excitation of LSPRs [22] 23 . As we predict here, 
the the metallic edge states on M 0 S 2 nanostructures can 
lead to the formation of highly localized ID plasmons - 
a result which adds edge-plasmonics to the list of this 
material’s unique properties with potential applications 
within nanoplasmonics and photocatalysis. 

In this communication, we use time-dependent density 
functional theory (DFT) to demonstrate the existence of 
a set of highly localized plasmons on the edges of an M 0 S 2 
nanoribbon, see Fig. [l] The fundamental properties of 
this new type of edge plasmon are investigated through 
a spectral analysis of the dielectric function which allows 
us to identify the plasmonic eigenmodes of the system in 
an unambiguous way. The deviations of the edge plasmon 
properties from those of an ideal ID plasmon are shown 
to arise from interband transitions between the metallic 
edge state and the M 0 S 2 conduction- and valence bands. 
The distinctly different plasmonic properties of the Mo- 
and S-edges are clearly seen in the spatially resolved loss 
spectrum suggesting that high resolution transmission 
electron microscopy could be used to probe the atomic 
structure of nanoribbon edges. 

All calculations were performed with the GPAW elec¬ 
tronic structure code, which is a grid-based implementa¬ 
tion of the projector-augmented wave method [24]. We 
considered a 1.5 nm wide sulfur terminated zigzag M 0 S 2 
nanoribbon, see Fig. [l] for a front view of the structure, 
and supplemental materials for a detailed image of the 
structure. This edge structure should be favorable at 
high chemical potential of sulfur, which can be tuned by 
the choice of sulfiding agent [25 . The ribbon was placed 
in a supercell including 20 and 12 A of vacuum in the two 
perpendicular directions. The wave functions were rep¬ 
resented in a plane-wave basis set with an energy cutoff 
of 340 eV. The structure was relaxed using the Perdew- 
Burke-Ernzerhof (PBE) exchange-correlation functional, 
while the GLLB-SC functional, including GLLB type ex- 
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Edge plasmons - Induced potential 



FIG. 1: Induced electric potential of the three plasmon 
eigenmodes at q = 0.2A , calculated as the 
eigenvectors of e q (r,r f ,u). The structure of the ribbon 
is shown in front view, with the ribbon running into the 
plane of the page. The Mo- and S-edge, covered by 
sulfur dimers, is found to the right and left of the 
ribbon respectively. 


change by Gritsenko et al. and PBEsol correlation [26] . 
was applied to calculate the single-particle states used as 
input for the linear response calculation. The GLLB-SC 
functional has been shown to improve the PBE band gaps 
of semi-conductors [26 and yield better descriptions of 
plasmon resonances in noble metals [27] due to its im¬ 
proved description of the d-band energies. However, for 
the M 0 S 2 nanoribbon GLLB-SC gives only minor correc¬ 
tions to the PBE band energies. The ID Brillouin zone 
was sampled with 128 k-points along the ribbon axis. A 
spin polarized calculation yielded zero magnetic moment 
for the structure. 

The dielectric matrix was calculated in the random 
phase approximation (RPA) following reference [28] : 

4tt 

e G ,G'(q,w) = <*G,G' - | q + Q| 2 XG,G / ( q ^)- (!) 

For the calculation of the non-interacting density 
response-function, x°> empty states were included up 
to 12 eV above the Fermi level, which was sufficient to 
converge the low energy plasmon peaks of interest. An 
energy-cutoff of 30 eV was used for the reciprocal lat¬ 
tice vectors G and Gb We used a non-linear frequency 
grid from 0 to 12 eV, with an initial grid-spacing of 0.01 
eV at uo = 0 and a smearing of 0.02 eV. The dense fre- 




FIG. 2: (a) Electronic bandstructure revealing three 
metallic edge-states (I-III). b) The three eigenvalues for 
the dielectric matrix fulfilling Ree n (g,ce) = 0 plotted for 
q = 0.06A . The color (red,blue,green) illustrates the 

one to one connection between the electronic bands and 
the plasmon eigenmodes. 

quency sampling and small broadening was necessary at 
low energies in order to resolve all intra-band transitions. 
The Wigner-Seitz truncated Coulomb approximation [29] 
was used in order to avoid interaction between supercells, 
which would otherwise give a significant energy shift at 
low momentum transfers. 

The calculated band structure of the nanoribbon 
shown in Fig. [ 2 ] reveals three metallic edge-states (labeled 
I-III) in agreement with previous calculations on a MoS 2 
nanoribbon with S-terminated edges [25]. The remaining 
bands produce a band-gap of approximately 2 eV. When 
the width of the ribbon is doubled to 3.3 nm the edge 
states are unaffected while the bandgap is decreased to 
1.8 eV which agrees with the value we find for the infinite 
single-layer MoS 2 . The Kohn-Sham wavefunctions I and 
III are almost completely localized on the S-dimers on the 
S- and Mo- edge of the ribbon respectively. The wave- 
functions are linear combinations of p^-orbitals, forming 
the one-dimensional metallic states along the edge. (Ill) 
has a large Fermi-velocity and is close to linear up to 1 eV 
away from the Fermi level. The wave function of state II 
resides at the outermost Mo-atom at the Mo-edge, and 
extends a few atoms into the ribbon. The number of 
edge states and their properties depends on the particu¬ 
lar edge-configuration. In this case the presence of sulfur 
dimers on the Mo-edge is crucial for obtaining the state 
(III). However, zigzag MoS 2 nanoribbons are in general 
found to be metallic opposed to armchair types [30 . We 
note that MoS 2 nanostructures could also have topolog¬ 
ically protected metallic edge states due to strong spin- 
orbit coupling in MoS 2 , which was not included in the 
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present study. In the following we show that the metallic 
states give rise to plasmon states at the ribbon edges. 

The plasmon modes of the system were investigated 
following a recently developed spectral analysis method 
for the dielectric function m- The eigenvalue equation 


X] e G,G'(q,w)</>£;,(q,w) = £n(q,u)<t>G.( q> w )> ( 2 ) 


is solved to obtain the frequency-dependent eigenvalues, 
e n , and eigenvectors, </> n , of the microscopic dielectric 
function. The plasmon energies are identified as those 
frequencies where the real part of an eigenvalue van¬ 
ishes: Ree n (q, cjp) = 0. When Ime n (q, cj) does not 
vary too much around the plasmon frequency, a vanish¬ 
ing real part coincide with a maximum in — Ime“ 1 (q, cj). 
The corresponding eigenvector </> n (q, ujp) gives the spa¬ 
tial form of the induced potential associated with the 
plasmon oscillation. 

For the M 0 S 2 nanoribbon we find three distinct plas¬ 
mon eigenmodes in the low-energy regime. The dielec¬ 
tric eigenvalue curves corresponding to these modes are 
for momentum transfer q = 0.2A 
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shown in Fig. 

Each mode arises from intraband transitions within one 
of the three metallic edge bands as indicated by the color 
code. The plasmons are clearly observed as peaks in 
Ime“ 1 (u;) shown in the lower panel. In contrast, the 
imaginary parts, Ime n , have peaks at the energy of the 
intra-band transitions associated with the edge-states 
and occur at qvp. The n = 3 plasmon has a higher en¬ 
ergy than the other two modes, partially due to the larger 
Fermi velocity of electrons in state (III). The induced po¬ 
tentials of the plasmon eigenmodes, given as the eigen¬ 
vectors of the dielectic matrix, are shown in Fig. [l] The 
edge plasmon corresponding to n = 3, is clearly local¬ 
ized on the S-dimers of the Mo-edge, and stems from the 
electronic state (III). The extend of the induced charge 
density of the edge plasmons are thus limited by the ex¬ 
tend of the electronic edge states, so that the plasmons 
will be atomically localized independently of the size of 
the M 0 S 2 ribbon. 

In Fig. [3a] the plasmon energies are plotted for in¬ 
creasing momentum transfers, and are clearly seen to 
be blue-shifted compared the the single-particle intra¬ 
band transitions (dashed lines) that disperse linear with 
q. In order to map out the effect of Landau damp¬ 
ing in (g, cj)-space, the sum of single-particle transitions, 
Ime n (g, cj) = ImTr e(q, cj), has been indicated in 
gray contours. The damping is seen become significant 
around 2 eV, where inter-band transitions between the 
valence and conduction band states of the M 0 S 2 sheet 
set in. 

At vanishing momentum transfers, the weight of single¬ 
particle transitions decreases significantly, since the num¬ 
ber of single-particle transitions scales as q. For an ideal 
ID electron gas the weight of the plasmon resonance 


should also go to zero in the limit q 0 [32]. How¬ 
ever, our calculations predict that the plasmon weight is 
almost constant in the considered g-range, and actually 
increases slightly for small momentum transfers. In or¬ 
der to explain this behavior, the result is compared to 
the pure intra-band plasmon, obtained by including only 
the single metallic band (III) (see Fig. [2| in the calcu¬ 
lation of XQ G ,(cj,q). This allows us to compare to the 
case of an undamped plasmon. As seen in Fig. |3b[ this 
approach yields a plasmon mode that is blue-shifted in 
energy compared to the full calculation and a weight that 
approaches zero for small q as expected for an ideal ID 
metal. Therefore, we conclude that the stagnant weight 
of the plasmon in the full calculation is due to damp¬ 
ing from coupling to interband transitions. The effect of 
interband transitions can be described by an effective di¬ 
electric constant of the medium that screens the Coulomb 
interaction : V(q) e^ er (g, u;)E(g). The screened in¬ 
teraction will be reduced in magnitude compared to V (g) 
and will have a finite imaginary part at energies corre¬ 
sponding to the interband transitions. These properties 
lead to a reduction of energy and life time of the full plas¬ 
mon compared to the bare intraband result. However, for 
very low momentum transfers, the plasmon is well sepa¬ 
rated from single-particle transitions, and a sharp reso¬ 
nance is obtained. Therefore the plasmon coincides with 
the undamped result, and is expected to have a long life¬ 
time in this regime. For the other two modes, we expect 
the damping by single-particle transitions to be lesssig- 
nificant due to separation in energy. 

Following reference [32]; a simple model for the plas¬ 
mon dispersion in the long wavelength limit can be ob¬ 
tained by approximating yo as that of a ID linear band 
with transitions at qvp: 


Xo(q,u) 


q \ 1 


7r 


uj — qvp + irj 


1 

uj + qvp + irj 


( 3 ) 


The dielectric function is obtained from the RPA expres¬ 
sion e = 1 — V(q)x ce where V(q) is the ID Coulomb 
potential. The inverse dielectric function can then be 
cast on the form 


e l (q,v) = 1 + 


a 

uj — c u(q) + ir] 


a 

uj + uj(q) + ir]' 


( 4 ) 


The pole of e 1 for positive frequencies and in the limit 
of 77 = 0 yields the plasmon dispersion 


uj(q) = qv F Jl + (5) 

V ttFf 

and the weight of the resonance equals a = 

gH(g)/^7ry / lT^^^. This result is identical to the 

eigen energies obtained from the interacting Tomonaga- 
Luttinger model with a ID Coulomb interaction. In 
this sense the RPA is exact for a ID system 133]. The 
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FIG. 3: (a)Energy dispersion of the three edge plasmons, showing plasmon energy (dots) as well as intra-band 
transitions (dashed lines). The dot size indicates the plasmon weight, defined as the integral over the loss peak. The 
gray contours map the regions of single-particle transitions, that contribute to damping of the plasmons. (b) 
Comparison of the dispersion of the S-dimer plasmon to the undamped intra-band mode obtained by only including 
the metallic state (III) in the response calculation. The inset shows the evolution of plasmon weight with q for the 
undamped and the full calculation. A ID RPA model valid at low q , identical to the Tamonaga-Luttinger result, has 

been fitted to the intraband results. 


Coulomb interaction for a thin wire of confinement a can 
be written for small q as V(q) = Vo|ln(ga)|, where Vo is 
a system dependent interaction strength accounting for 
e.g. the screening due to interband transitions. The ex¬ 
pressions for the plasmon frequency and spectral weight 
derived from the linear band RPA model were fitted to 
the pure intraband results in Fig. [3bJ As expected from 
the highly linear form of band (III), the quality of the fit 
is excellent. From the fit we obtained a = 1.34 A which 
is in the same order as the radius of the orbitals of the 
edge state. Comparing the dispersion and weight of the 
model to the results of the full ab initio calculation, we 
find agreement for small (g,a;), after which the plasmon 
energy and, in particular, the plasmon weight start to de¬ 
viate in a nonlinear way. This shows that the screening 
varies as a function of q and uj as is generally the case, 
such that the interaction strength cannot be written as 
a constant Vo as assumed in the model. 

Experimentally, plasmon excitations can be probed by 
electron energy loss spectroscopy (EELS). When per¬ 
formed with a highly confined electron beam, the spa¬ 
tial form of the loss spectrum can be obtained with few 
Angstrom resolution [34] . 

The loss spectrum is defined as the power dissipated 
in the medium due to an external potential, (j) ext (r)e' lujt : 

PH-/ J drdr'(t> ex t(r)x2(r,r , ,uj)(t> ex t(r'). (6) 

Here \2 is the imaginary part of the interacting density 


response function. In the case of EELS, the external 
potential corresponds to the Coulomb potential of a fast 
electron moving at constant velocity v e emitted at point 

1*0: 


( fiext (C t) 


47re 2 

|r - r 0 — v e t| ’ 


(7) 


The final expression for the loss spectrum is given in the 
supplemental material [35,. In Fig. [4] the spatially re¬ 
solved EELS spectrum is plotted for two different values 
of q x for a beam position that is varied acroos the ribbon. 
The spectrum is dominated by the S-dimer plasmon on 
the Mo-edge, which reflects the large coupling strength 
of this mode. The other two edge plasmons are also 
visible in the spectrum, particularly at larger momen¬ 
tum transfers where the Coulomb potential is less long 
ranged, which results in more localized features. Further¬ 
more, the modes can be distinguished by their difference 
in energy and dispersion with q. We suggest this ap¬ 
proach could be applied to verify the existence of the 
edge-plasmons in M 0 S 2 nanoribbons or clusters. 

In conclusion, our first-principles calculations predict 
the existence of a unique type of plasmon localized on 
the edge of a zigzag M 0 S 2 nanoribbon. The fundamen¬ 
tal properties were investigated in detail and its clear 
signature in the spatially resolved electron energy loss 
spectrum was demonstrated. Finally, we proposed that 
the edge plasmons could be utilized to tune the photo- 
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FIG. 4: Spatially resolved electron energy loss spectrum 
calculated for two different values of momentum 
transfer, q x . We observe a strong signal associated with 
the S-dimer plasmon at the Mo-edge, while the two 
other plasmon modes give a weaker response. At low 
momentum transfers the signal is less localized due to 
the long range of the coulomb potential. 


catalytic activity of M 0 S 2 nanoparticles as was recently 
demonstrated for metal nanoparticles {22,, 23 . 
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